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Strongly correlated phases exhibit collective carrier dynamics that if properly harnessed can enable novel 
functionalities and applications. In this article, we investigate the phenomenon of electrical oscillations in a 
prototypical MIT system, vanadium dioxide (VO2). We show that the key to such oscillatory behaviour is the 
ability to induce and stabilize a non-hysteretic and spontaneously reversible phase transition using a 
negative feedback mechanism. Further, we investigate the synchronization and coupling dynamics of such 
VO2 based relaxation oscillators and show, via experiment and simulation, that this coupled oscillator 
system exhibits rich non-linear dynamics including charge oscillations that are synchronized in both 
frequency and phase. Our approach of harnessing a non-hysteretic reversible phase transition region is 
applicable to other correlated systems exhibiting metal-insulator transitions and can be a potential 
candidate for oscillator based non-Boolean computing. 



Rhythmic activity is ubiquitous in the physical world; and the process of synchronization of coupled 
oscillatory systems can be observed in the circadian rhythms of organisms^"^, the synchronization of an 
ensemble of fireflies^ and is a distinguishing feature in certain structures of our brain including the 
neocortex and the thalamus^ Participating elements of a coupled system have the fundamental ability to exhibit 
a non-linear response to forcing stimuli; and hence under certain conditions, a collectively ordered state emerges 
from the complex dynamics of such a non-linear system. This motivates the development of physical systems of 
synchronous oscillators which are energy efficient, compact, scalable and could potentially enable large scale 
biologically inspired computing algorithms such as associative computing for applications like real-time image 
and pattern recognition^"^^. 

In this paper, we investigate the dynamics of self-sustained electrical oscillations in a prototypical correlated 
system vanadium oxide (VO2) which undergoes an abrupt first order metal-insulator transition with up to five 
orders of change in conductivity. Further, we explore the possibility of realizing a coupled system of such VO2 
based relaxation oscillators by exploring their synchronization and coupling dynamics. Through experiment and 
circuit simulation, we demonstrate that a system of capacitively coupled VO2 oscillators exhibits charge oscilla- 
tions that are synchronized in frequency and phase (anti-lock), thereby allowing the opportunity to harness the 
non-linear dynamics associated with their synchronization in computational algorithms. Our work enables a 
unique application and design space for a system of coupled correlated oxide based oscillators particularly suitable 
for applications in the domain of non-Boolean computing^. 

Results 

Single Oscillator Dynamics. Electrical oscillations in VO2 harness the first order metal-insulator phase 
transition, which in this case is triggered electrically. It is well accepted that these oscillations represent a 
triggering of the insulator- to -metal transition (IMT) followed by a resetting metal-to-insulator transition 
(MIT) leading to spontaneous oscillations^l However, the physical origin of the electrically driven phase 
transition and therefore, the fundamental origin of the electrical oscillations in VO2 is widely debated as being 
electronic^^"^^ or electro -thermally driven^°'^\ Researchers have also shown evidence of electric field induced 
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Figure 1 | Stabilizing a non-hysteretic and reversible phase transition in VO2 and single oscillator dynamics, (a) Schematic of the circuit used to 
access the non-hysteretic transition regime in VO2 using the negative feedback generated by the series resistor Rs. (b) Rs modifies the VO2 phase transition 
dynamics through negative feedback. A non-hysteretic reversible (NHR) transition regime can be accessed using the appropriate Rs (>Rc). Electrical load 
lines for non-hysteretic ( 1 ) and hysteretic (2) switching are superimposed, (c) Time domain waveform of the VO2 relaxation oscillator. (Inset) shows the 
waveform for one oscillation period with the two time constants (xi = 0.35 [is and X2 = 0.65 |as) associated with IMT and MIT. (d) Power (Mean Square 
Amplitude in dB) spectrum of a single oscillator. The fundamental frequency of the oscillator here is ~1 MHz. (e) Simultaneous measurement of the 
voltage and scattered x-ray diffraction intensity for an oscillating VO2 device. Based upon the Bragg peaks, the structural phase of VO2 is labeled. (Inset) 
shows the [40 2] and [002] Bragg peaks of the Ml and rutile phase ofV02, respectively. The dashed line indicates the 0/20 angle ( = 51.714°) at which this 
measurement was performed. Details of x-ray diffraction measurement are described in the Supplementary S4. (f) Frequency scaling of the VO2 relaxation 
oscillator with R^. The frequency increases as the external Rs is scaled down. The device resistance has also been scaled to further increase the output 
frequency. 



nucleation^^ as the origin of IMT in VO2, similar to switching in 
chalcogenide glasses^^'^^. Some reports elucidate the formation of 
conducting filaments^°'^^ that trigger the IMT. This debate is linked 
to the larger controversy over the fundamental driving physics of the 
phase transition in VO2, particularly whether it is a Mott^^'^^'^^'^^ or a 
Peierls^^"^° transition or a weighted contribution from both the 
mechanisms^ \ Despite the uncertainty over the unequivocal origin 
of the electrically driven transition in VO2, our work, as well as other 
reports in literature^^'^^ show that the dynamics of such relaxation 
oscillators can be modulated with temperature, total circuit 
resistance and capacitance. 

To understand the single oscillator dynamics, we first explore the 
direct current (DC) characteristics of the oscillator circuit which 
consists of a two terminal VO2 device (see Methods) with a series 
resistor Rs (Fig. la). The dimensions of the device used for measure- 
ments (except in experiments relating to Fig 2b, c) are Lvo2 — 4 jim; 
Wvo2 = 10 jim. The device size used in experiments relating to 
Fig. lc,d is Lvo2 = 4 jim; Wvo2 = 80 jam. The series resistor pro- 
vides a negative feedback and, therefore, modifies the transition 
dynamics. Unless stated otherwise, all electrical measurements 
reported are performed near room temperature (18°C). 

When electrically driven across the phase transition with Rs = 0 Q 
(Fig. lb), the IMT and the MIT in VO2 occur at critical fields, E2 and 
El, respectively. A current compliance is set to limit the metallic state 
current and prevent permanent device breakdown. While the exact 
magnitude of these electric fields depends on the origin of the phase 
transition (thermal/electronic), they represent a quantitative mea- 
sure of the stimulus threshold that enables the IMT and the MIT to 
occur. If the transition is thermally driven then E2 will represent the 



electric field at which there is enough current density and corres- 
ponding thermal flux (most likely to be highly localized) to induce 
the IMT. Similar consideration applies to electronic case as well as to 
El. Therefore, irrespective of the specific physical origin of the trans- 
ition, the state of VO2 can be described by knowing the electric field 
across the device. E2 represents the field at which the conducting 
metallic phase can be stabilized to form a low resistance state while 
El is the field at which the metallic phase can no longer be stabilized 
and the device returns to its high resistance insulating state. Even 
though the fundamental mechanism of electric field induced switch- 
ing remains nontrivial, the delay time^^ as a function of DC trigger 
voltage as well as the temperature dependence of these critical field 
points is shown in the Supplementary SI. These critical fields differ in 
magnitude giving rise to a hysteresis in the electrically driven phase 
transition. 

In contrast to this hysteretic transition, a non-hysteretic transition 
could be induced by modifying the phase transition dynamics such 
that, triggering the IMT (by surpassing E2) results in a negative 
differential resistance (NDR) characterized by an increase in con- 
ductivity and a simultaneous reduction in the field across the VO2 
device. If the electric field across the VO2 now drops below Ei, it 
makes the metallic phase unstable resulting in a spontaneous MIT. 
Representing this scenario through an electrical load line (Fig. lb), a 
VO2 device operating on load line '1' will undergo a non-hysteretic 
transition, whereas a device operating on load line '2' will go through 
a hysteretic transition. 

While the NDR and the load lines described above cannot be 
achieved using the Rs = 0 Q configuration, the addition of a simple 
Rs in series makes this feasible. Rs provides a stabilizing negative 
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feedback (see Supplementary S2) compensating the positive feedback 
associated with the abrupt change in conductivity. This negative 
feedback modulates the current- electric field characteristics across 
the VO2 and establishes the load line relation: 



Evo2 — E2 — Alvo2Rs/Lvo2 



(1) 



where Alvo2 — Ivo2(Evo2) ~ Ivo2(E2) and Lyo2 is the VO2 channel 
length (see Supplementary S3 for derivation of equation (1) and (2)). 
From equation (1), it is evident that Rs tunes the VO2 current- electric 
field characteristics and sets the criteria for non-hysteretic/hysteretic 
switching. The critical resistance, Rc, that enables the non-hysteretic 
regime of operation, ensures that the metallic phase is never com- 
pletely stabilized and is given by: 



Rc> 



Wv02tv02 



-El 



(2) 



where Wyo2 (^^(^ W02 are the VO2 channel width and thickness, 
respectively, and a^i and a^^ are the equivalent insulating and metal- 
lic state conductivities, respectively. 

The conductivity in the NDR regime as indicated by Al yo2 can be 
tuned by the electric field (equation (1)) with conductivity increasing 
as the electric field drops (and vice-versa). Such a system can oscillate 
if the conductivity and electric field across the device modulate each 
other in a way that the restoring force enables the system to trace the 
same electrical trajectory periodically^^. This criterion is fulfilled only 
in the non-hysteretic operating regime and hence oscillations are 
observed only when Rs > Rc- 

The time domain waveform of the VO2 oscillator (Fig. Ic; single 
time period shown as inset) shows an exponential voltage buildup 
and decay (across Rs; 5 kQ in this case) associated with MIT and the 
IMT in VO2 respectively and is characterized by two R-C time con- 
stants Ti and T2 (see inset Fig Ic) which primarily control the oscil- 
lation frequency. Xi ( = 0.35 jis), associated with the IMT is smaller 
than T2 ( = 0.65 |is). These time constants are much larger than the 
intrinsic electronic switching time^^ of VO2 and therefore oscillation 
frequency can be scaled by reducing R and C through VO2 device 
scaling. The corresponding power spectrum for the oscillations is 
shown in Fig. Id. The exponential nature of the time domain wave- 
form results in a gradual decay of the frequency harmonics and 
therefore the first harmonic is approximately 10 dB below the fun- 
damental frequency. 

In-situ nano X-ray diffraction (see Supplementary S4) performed 
on the VO2 film stabilized in an oscillating state (Fig. le) confirmed 
that a structural transition occurs during each oscillation cycle sim- 
ilar to the DC electrically driven phase transition^^. The low resist- 
ance regime of VO2 during the oscillation cycle is characterized by 
the rutile metallic phase, whereas the insulating regime is character- 
ized by the monoclinic Ml phase. This is expected for the —0.9% 
compressively strained VO2 films at room temperature^^ (see 
Methods). 

We then study the oscillator frequency (fosc) and dimensional 
scalability which is crucial to implementing dense, large scale VO2 
oscillator circuits. Since fosc is set by the R-C time constant of the 
oscillator circuit, we explore fosc scaling with Rs (Fig. If). As Rs is 
reduced, frequencies up to 1 MHz are experimentally realized. All 
experiments in this work are performed using discrete components 
(resistors) and eliminating the parasitic elements associated with the 
present wiring schemes through monolithic integration along with 
device scaling will further increase fosc as observed with other corre- 
lated materials^^'^^. 

Further, we note that Nardone et aV^ who demonstrated relaxa- 
tion oscillators in chalcogenide phase change memory, have illu- 
strated some of the stochastic features associated with such 
oscillations including the observation of multiple oscillatory patterns 
such as oscillations with decaying amplitudes. We have experiment- 
ally observed sustained stable oscillations over 2.5 X 10^ cycles (see 



supplementary SI) with no decaying amplitudes. Additionally, it is 
expected that the system of coupled VO2 oscillators (described 
below) will have higher stability to stochastic variations due to 
mutual feedback induced stabilization^^. 

Coupled VO2 oscillators and their dynamics. The realization of 
inductor free and scalable non-linear oscillatory systems using 
VO2 can be better exploited in the application domain if the 
oscillators can be synchronously coupled giving access to the large 
array of non-linear dynamics which are associated with their 
synchronization. Additionally, long range electrical charge based 
coupling is crucial to implementing large oscillator arrays com- 
pared to other state- variables (like spin), which may have limited 
coherence length at room temperature. 

To investigate the coupling dynamics, we first explore 'injection 
locking' wherein a single oscillator is first excited with an external AC 
sinusoidal signal coupled through a coupling element which in this 
case is a capacitor Cc (Fig. 2a). The oscillator locks onto this external 
frequency only when it is close to the oscillator's natural frequency 
and within its lock- in range so that a large power output is observed 
at the input signal frequency (Fig. 2b). The ability of a single oscillator 
to lock onto an external AC signal implies that the two coupled 
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Figure 2 | Injection locking of VO2 oscillators, (a) Schematic diagram of 
the injection locking circuit wherein an external AC signal is coupled to the 
oscillator through a coupling capacitor Cc- (b) Power output 
characteristics of a VO2 relaxation oscillator (Rs = 68 kQ; —lORc) 
capacitively coupled to an external frequency signal (Cc =1.15 nF) 
showing locking behavior for frequencies close to its natural oscillation 
frequency. 
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Figure 3 | Synchronized oscillations in VO2 oscillators and their coupling dynamics, (a) Schematic of two capacitively coupled VO2 oscillators. 
Individually, the oscillators are tuned to different oscillating frequencies using Rgi and Rs2 of 38 kQ ( = 6Rc) and 47 kQ ( = 7.5Rc), respectively, 
(b) Experimental and simulated power spectrum of both oscillators before and after coupling. The individual oscillator frequencies lock to a new 
frequency when coupled capacitively (Cc = 680 pF). (c) Effect of Cc on the full width at half maximum (FWHM) of the spectral line width and resonant 
frequency of coupled VO2 oscillators. Decreasing FWHM with increased Cc reflects the increased coupUng strength, (d) Effect of external resistance 
symmetry on the coupling strength. The coupling strength is maximal when the oscillators are identical (equal Rg). 



oscillators could mutually lock and exhibit synchronized oscillations 
as investigated further. 

Fig. 3a shows a schematic of the coupled oscillator circuit. Prior to 
coupling, the two individual VO2 oscillators have fosc of 48.43 kHz 
and 37.27 kHz, respectively (Fig. 3b). fosc is set by choosing appro- 
priate Rs values: 38 kQ ( = 6Rc) and 47 kQ ( = 7.5Rc), respectively. 
After the oscillators are capacitively coupled (Cc = 680 pF), they 
frequency lock and converge to a single resonant frequency 
(22.95 kHz). This is accompanied by a sharp narrowing of the spec- 
tral line width (Fig. 3b). The reduced line- width after locking reflects 
the mutual feedback induced stabilization against noise^^, implying 
that the synchronized state is a stable configuration for the coupled 
system. 

The use of capacitive coupling is motivated by the high pass filter- 
ing characteristics of the coupling configuration, which ensures that 
the individual oscillators can synchronize without mutually affecting 
their DC quiescent point. Cc modifies the frequency dynamics of 
each oscillator through an effective coupling term Cgff 

Ceff = CcCi+CcC2 + CiC2 (3) 

where Ci and C2 are the net capacitances of each oscillator (see 
macro-model in supplementary). The Cc element stores and dynam- 
ically redistributes reactive (non-dissipative) power between the 
oscillators facilitating synchronization. Further, Cc also tunes the 
oscillator coupling dynamics and influences the strength of the coup- 



ling interaction between them as indicated by the spectral line-width 
(full width at half maximum; FWHM) (Fig. 3c). Additionally, R^, 
which tunes the natural frequencies of the individual oscillators, also 
affects the strength of the coupling between the two oscillators. In the 
symmetric oscillator configuration (equal Rs, identical frequencies 
for both oscillators), the strongest coupling is observed (Fig. 3d). The 
ability to tune the synchronization dynamics of the oscillators is 
crucial from an application standpoint. The experimentally observed 
frequency locking is also verified through circuit simulations using 
macro -models (see Supplementary S5), which describe the system's 
time evolution (Fig. 3b). 

The synchronization dynamics of this coupled oscillatory system 
are analysed using the equivalent model shown in Fig. 4a. The inter- 
action between the two oscillators, each of which traverses between 
the metallic (M) and the insulating (I) phases, is analysed through the 
resulting voltage flow diagrams. Figure 4b illustrates the flow dia- 
grams of the coupled systems as they evolve over time. As the VO2 
device in each oscillator can be either metallic (M) or insulating (I), 
the resultant phase plane for the two VO2 devices will consist of four 
quadrants where the oscillators can be in state MM, MI, IM or II. 
After coupling, the system will settle on one fixed orbit/trajectory, 
which may go through these four quadrants. As the device and circuit 
parameters are changed, the stable orbit will change indicating a 
different path in the phase plot. The blue arrows indicate a family 
of such orbits corresponding to the time evolution of the system 
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Figure 4 | Analyzing the coupling dynamics of VO2 oscillators, (a) Equivalent circuit model for capacitively coupled VO2 oscillators. Ci and C2 
are the effective capacitances of each oscillator, which include the internal VO2 capacitance associated with VO2 film and the device (CV02) as well as 
parasitic capacitances. Rp 00. The same circuit is used to simulate the experimentally observed frequency locking and has been verified with 
experimental data as shown earlier, (b) Voltage flow diagram in phase space for steady state coupled oscillations. 1-2-3-4 is the electrical trajectory of the 
capacitively coupled configuration. It reveals that the coupled oscillators primarily move in the M-I and I-M phase planes, (c) The electrical trajectory of 
the coupled oscillator configuration in phase space, (d) Simulated time domain waveform of the two VO2 oscillators after coupling (top). The to and fro 
movement of charge (Qc) across the capacitor enables the two oscillators to synchronize (bottom). 



shown in Fig. 4a. The stable orbit corresponding to our experiment 
has been shown in red. The stable orbit (red curve in Fig. 4b) of the 
system has been exclusively shown again in Fig. 4c as a phase plot (Vi 
vs. V2) along with the corresponding metallic and insulating phases 
for each VO2 device (the numbers 1,2,3,4 show the related arm in the 
flow diagrams in Fig. 4b). The phase diagram in Fig. 4c reveals that 
the synchronized oscillators primarily move between the M-I and I- 
M planes such that as one VO2 device transitions into the metallic 
phase, the other device evolves into the insulating phase (b ^ c; d ^ 
a). This is also marked by charging/discharging of the coupling capa- 
citor Cc- The region a ^ b; c ^ d defines the short phase space 
trajectory where both the devices have the same phase and represents 
the time during the coupling capacitor Cc reverses polarity. The 
corresponding time domain waveform of the two oscillators after 
coupling is shown in Fig. 4d. As the oscillators move in the phase 
space, the periodic sourcing/draining of charge (Fig. 4d) by the coup- 
ling capacitor Cc enables the oscillators to synchronize. 

Discussion 

The non-hysteretic transition dynamics that enable such spontan- 
eous electrical oscillations in VO2 inspires the exploration of other 
correlated systems for similar applications^""^^. The coupled oscil- 
latory dynamics demonstrated here provides an important mile- 
stone of experimental demonstration of coupled synchronous 
oscillations in correlated oxide materials and unravels new oppor- 



tunities towards exploring large coupled oscillatory systems. We 
anticipate a wide array of dynamics in larger arrays of such coupled 
oscillators that could potentially be harnessed to develop a robust 
hardware platform capable of supporting novel computing schemes 
and algorithms. 

Methods 

The VO2 two terminal devices are fabricated on a 17.7 nm thick VO2 films. The films 
are epitaxially grown on a Ti02 (001) substrate employing reactive oxide molecular- 
beam epitaxy using a Veeco GEN 10 system. Due to epitaxial mismatch, the VO2 films 
are biaxially strained in compression by —0.9%. First the electrodes are patterned 
using contact lithography followed by electron beam evaporation of Pt/Au (15 nm/ 
80 nm) and lift-off in RemoverPG at 70°C. Next the channel width and isolation are 
defined by electron beam lithography followed by a CF4 dry etch. Finally, the resist is 
stripped with RemoverPG at 70°C. 
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